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Abstract 

In this letter, a coupled system of viscous Burgers' equations with 
zero Dirichlet boundary conditions and appropriate initial data is con- 
sidered. For the well-known single viscous Burgers' equation with 
zero Dirichlet boundary conditions, the zero equilibrium is the unique 
global exponential point attractor. A similar property is shown for the 
coupled Burgers' equations, i.e., trajectories starting with initial data 
which is not too large approach the zero equilibrium as time goes to 
infinity. This "approaching" or convergence is not necessarily expo- 
nentially fast, unlike the single viscous Burgers' equation. 
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1 Introduction 



We consider the following coupled Burgers' equations 

ut = Uxx - uux - a{uv)x, 



vt = Vxx - vvx - b{uv)x, 



together with Dirichlet boundary conditions 

u{0,t) = u{l,t) = 0, 
v{0,t) = v{l,t) = 0, 

and appropriate initial conditions 

u{x,0) = uoix), 
v{x,0) = vo{x), 



(1.2) 



(1.3) 



for 0<x<l,t>0. Here a, b are constants. 

This coupled system, derived by Esipov , is a simple model of sedimen- 
tation or evolution of scaled volume concentrations of two kinds of particles 
in fluid suspensions or colloids, under the effect of gravity. The constants 
a, b depend on the system parameters such as the Peclet number, the Stokes 
velocity of particles due to gravity, and the Brownian diffusivity. In ^ Es- 
ipov reported numerical simulations for (|1.1|)- (|1.3| ) and compared the results 
with experimental data. 

In this letter, we consider dynamical aspect of the coupled system. We 
show that trajectories (orbits of solutions) of this coupled system approach 
the zero equilibrium as t ^ +oo when the initial data u{x, 0), v{x, 0), is not 
too large in some Sobolev norm. That is, the w-limit set of these trajectories 
is the zero equilibrium. 



2 Single viscous Burgers' equation 

We first recall an interesting property about w-limit set for the well-known 
(single) viscous Burgers' equation (||2|), i.e. equations in (LI) without the 
nonlinear coupling term 

Wt = Wxx - WWx, (2.1) 
with boundary and initial conditions 

w{0,t) = wil,t) = 0, 
w{x, 0) = wo{x). 
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The interest in the Burgers' equation arises because it is a simple one dimen- 
sional analog of the Navier-Stokes equation. The importance of the Burgers' 
equation is due to the nonlinear convection term uux- 

In the following, L'^{0, 1), L~(0, 1), H^{0, 1) and H^{0, 1) are the usual 
Sobolev spaces, while C(0, 1) is the space of continuous functions. We denote 
by II • II the usual -L^(0, 1) norm. All integrals / are with respect to x G [0, 1], 
unless specified otherwise. 

Multiplying the equation ( |2.1D by w and integrating over x £ [0, 1], we 

get 

II l|2 II ||2 ^ II ||2 /r, o\ 

^J^M =-\M <-\\w\\ , (2.3) 
where we have used the Poincare inequality 

/ w'^dx < f wldx, (2.4) 
Jo Jo 

in the last step. Thus by the Gronwall inequality (Q, p. 88), we further have 

\\w{x,t)f <\\w{x,0)\\'^e-'^\ t>0. (2.5) 

This means that all trajectories converge in the L^(0, 1)- norm to the zero 
equilibrium exponentially fast, i.e., the w-limit set (Q) of every trajectory 
is the zero equilibrium. The zero equilibrium is the global point attractor. 
See and for further results in this regard. This property also holds 
for higher dimensional viscous Burgers type convection-diffusion equations 
(§)• 



3 cj-limit set of the coupled viscous Burgers' equa- 
tions 

In this section, we show that trajectories of the coupled viscous Burg- 
ers' equations ( |1.1| )- (|1.3| ), whose initial data is not too large in Hq{0,1)- 
norm, converge to the zero equilibrium in the max-norm. This convergence, 
though, is not necessarily exponentially fast, unlike the situation for the 
single viscous Burgers' equation. 
Let 

m = \{\\u\\' + \\v\\% (3.1) 

J(t) = \i\\ux\\^ + \\vx\\^), (3.2) 
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Kit) = ^i\\u,,\\^ + \\v,,\\^). (3.3) 

For initial data u{x,0),v{x,0) G Hq{0,1), local (-in-time) existence and 
uniqueness of classical solutions to ( |l.lD -(pr^) can be easily shown by the 
usual semigroup method; cf. Theorems 3.3.3, 3.3.4 and 3.5.2. Global 
existence and uniqueness of classical solutions then follows once we show 
a priori that the solutions do not become unbounded in Hq{0, l)-norm at 
finite positive time. Moreover, the solutions, whenever they exist, are in 
H^{0, 1), /c = 1, 2, • • •; (cf. |], p.73). Note that I{t) < J{t) < K{t) whenever 
u,v exist, by the Poincare inequality. 

We will show that J{t) — > as t — > 0, when J(0) is bounded by some 



constants depending only on the system parameters a, 6 in (1.1). 

Multiplying the first and second equation in ( |1.1| ) by —Uxx and —Vxx, 
respectively, adding the two resulting equations and then integrating over 
X £ [0, 1], we obtain 



dJ 
Itt 



-ulx + uUxUxx + a{uv)xUxx - vlx + vVxVxx + b{uv) 



-ulx - Y - f^L - y + (a + b){uv)x{Uxx + Vxx) 
< / -^2 _ ^ _ 2 _ ^ 



+ (|a| + |6|)(y {UxV + UVxf)^{j {Uxx + Vxx?)^ 

< -2K-J^-J^f + {\a\ + mj2{ny + n'vl))H2K)'^ 



< -2K 




+2{\a\ + \b\)U ui v^ + J (3.4) 
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Now we further estimate the right hand side of ( |3.4D term by term. 
From the fact that 

= 2 I uurdx < 2\\u\ \ ■ \ \u^ 







for u £ Hq{0, 1), we get the so-called Agmon inequality 

< 2||n|| • llii^ll, (3.5) 
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where ||m||oo is the L°°-norm. 

Observe, using the Cauchy-Schwarz inequahty (1^], p. 183) and the above 
Agmon inequahty 

u'tdx = / uidu = —2 / uUxUxxdx 



1 \ i 



< 2\\u\\oo{j ul)'2{j ulS^ 

< 2V2{j u')Hj uDHj uDHj 

= 2V2{ [ u')H [ ul)H [ ulJ 



13 1 
< 8/4J4K2 

. 1 



< 8JK2 (3.6) 

In the last step, we have used the fact that I{t) < J{t). The same inequahty 
holds for / v^dx. 



Similarly, 



uj.dx = / u^du = —3 / uu^Uxxdx 



< 3\\u\\^ij ut)2{J ui^)-2 

< 3V2{ [ u')H f ul)H f ut)H I ui^r^ . (3.7) 



Thus 



{Jutdxf^ < 3V2{J u')kj uDkJ ulJ 

< 6V2J^K^2. (3.8) 



The same inequality holds for / v^dx. 
Moreover, note that 



rx 

= { I 2uUxdx)'^ — 4 / u'^(ix I u^dx. 



So we have 



J u^dx <4j u^dx j uldx < 16IJ < 16 (3.9) 
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The same estimate holds for / v'^dx. 



Substituting ( p^ ) , (|3.8D and into (|3.4|) , and using the Young 's in- 



equality {M, p. 108), we finally get 



dJ 



, < -2K + 8JK^ + 8V32^(\a\ + \b\)jiK'i 
dt 

< -2K +'-K + 164J2 + + 4^[8^/32i (|a| + \b\)]^j\?,.m) 
for e > 0. Taking e = | and noting that J(t) < K{t), we now have 

^ < -K + 40J^ + 36000(|a| + 
at 

< -J + 40J2 + 36000(|a| + |6|)^J^ (3.11) 
The comparison equation 

^ = - J + 40J2 + 36000(|a| + \b\)^J^ = f{J) (3.12) 
at 

has fixed points and 



_ -1 + VI + 9000( a + b 
J+ — — 



1800(|a| + |6|)4 

The third fixed point is negative and is discarded since J{t) is always non- 
negative by definition (|3^) . We calculate /'(O) < and /'(J+) > 0. So fixed 
point is stable while J+ is unstable; see ||l^, p. 187 or [|ll|, p. 8. That is, for 
the comparison equation ( p.l2| ), J{t) — > if J(0) < J+, while J{t) — > oo if 
J(0) > J4.. Due to the standard comparison result for ordinary differential 
inequalities and equations (|12|, P.69), the J{t) satisfying the differential 



inequality ( ^.11 ) also approach zero as t goes to infinity, when J(0) < J4 

We have thus shown that, if the initial data satisfies ^(||^io(-^)lP"'"ll''^o(-^)lP) < 
J+, the corresponding classical solutions and hence trajectories exist for all 
t > 0, since the i?o(0,l) -norm, i.e., J(t), in this case is bounded. More- 
over these trajectories approach the zero equilibrium as t ^ co. The zero 
equilibrium is the w-limit set of these trajectories. We remark that this 
convergence is not necessarily exponentially fast, unlike the single Burgers' 
equation. If |(||uo(x)|p + ||fo(a^)|P) > however, we cannot conclude any 
thing about the corresponding trajectories, based on the above dynamical 
system style analysis. 
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Note that \\ux\\, \\vx\\ are actually ^^^o(0, l)-norms of u and v due to the 
Poincare inequality, and note also that Hq{0, 1) is embeded in C(0, 1). So 
the above convergence of trajectories holds in the max- norm in C(0, 1). 

Therefore we obtain the following theorem 

Theorem 1 Assume that uq{x),vq{x) are in Hq{0,1) and satisfy 

2(ll%(^)ll +11^0(^)11 )< 225(|a| + |6|)4 
Then the global unique classical solutions of the coupled system ^^T^-l^^- 



(1.5) exist, and the corresponding trajectories (orbits of solutions) approach 
the zero equilibrium in max-norm. That is, the zero equilibrium is the oj- 
limit set of these trajectories. 
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